Compact binary waveform recovery from the 
cross-correlated data of two detectors by matched 
^ ; filtering with spinning templates 

j>^; L. Vereb^'2, Z. Keresztes^'^*, P. Raffai^^ Zs. Udvari\ M. Tapai^^ 

C^ ■ L. A. Gergelyi^t 

Si 

rvj . Abstract. We investigate whether the recovery chances of highly spinning waveforms by 

^-H ' matched filtering with randomly chosen spinning waveforms generated with the LAL package 

are improved by a cross-correlation of the simulated output of the LI and HI LIGO detectors. 
' ^' . We find that a properly defined correlated overlap improves the mass estimates and enhaces the 

^^' recovery of spin angles. 
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O ! 1- Method 

'""' ■ Recovering a gravitational wave pattern in the noisy detector output is a difficult problem. In 

^ , this work we present our attempt to recover a gravitational spinning waveform hinj immersed 

'kj^ \ in S5 LIGO-like noise by a modified matched filtering method. We calculate the overlap |lj 

^ • between the noisy injection hn,i and spinning templates htempiate as 

a ! 

Or, 7 1 \i^n,i\ri'templatel /-.\ 

["-n.i) i^template] — / ,, , , , ., , , ^ ■> \'-j 

Y \'>'n,i\'T'n,i) \i^template\'^templatel 

with _ _ 

where tilde and star denote Fourier transform and complex conjugate, respectively, and Sn (/) 
is the power spectral density of the noise. We choose /min = 50Hz and /max = 600Hz, lying 
in the best sensitivity band of the LIGO detectors, also the post-Newtonian prediction for the 
waveform is quite accurate there. 

On theoretical grounds, infinite long data series would be required for exact determination of 
the power spectrum. In order to achieve stability of the power spectral density of order of 1%, 
we would like to have at least 100 periods of the lowest frequency, therefore a minimal length of 
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Figure 1. The polarization /i-|- (left) and hx (right) of the spinning waveforms. 



the templates of 2 sec was imposed. We compute the above overlap both for the Hanford and 
Livingston detectors {O^ and O^, respectively). 

We also defined a correlated match, the overlap O^^^^, using the correlations of the Hanford 
and Livingston signals for /i„^j, htempiate and noises. (The correlation of hi and /12 in Fourier 
space is defined as H (/) = hi (/) /i2 (/)•) 
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Figure 2. Spinning waveforms at Hanford (top left) and Livingston (top right), mixed in LIGO 
S5-like noise (bottom). 
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Figure 3. Characteristic unfiltered (left) and filtered (right) LIGO noises. 



Table 1. The parameters of the gravitational wave polarizations (/i+, hx) of the injected signal. 
Masses mi, magnitude of the dimensionless spins Xii initial direction of the spin vectors given 
by cosKj and ■01 in the frame with the line of sight on the z-axis, and the initial orbital angular 
momentum Ln in the x-z plane, spanning the angle Q, and the distance (Il of the source. The 
initial phase is 0; the initial time being also fixed to 0. 



name 



mi(M(7 



Q) 



"12 (Mo) 



Xi 



X2 



COSKi 



^P 



COSK2 



■02 



e 



di (Mpc) 



injection 3.553 



3.358 0.983 0.902 0.984 1.109 0.978 0.957 1.430 



1 



Table 2. The angles 9, (p and ip (polarization angle) give the relation between a fictitious 
Earth-centered detector and the source frame. For the actual detector positions further rotation 
have to be taken into account [5]. These angles are also necessary for computing the antenna 
functions (F+,Fx). The gravitational signal is /i = /i+i*+ + fixF^. 
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injection 


3.657 


0.278 


0.000 



For the two detectors we also employed the quantity O^^, according to Ref. [2]: 

\'T'n,i\'T'template) j{ ~r \iT"n,'i\i^templatelL 



0^^ = 



a/ \\iT"n,i\iT'n,i/ l{ ' \'^n,i\'^n,i) l) \\'^tem,plate\'^templatel u + {'T'templatel'T'template) l) 



(3) 



Therefore we could compare four type of matches. 

The spinning waveforms, based on Ref. [3], were computed by using SpinTaylor Ref. [4j, while 
the antenna functions based on Ref. [S] with the XLALComputeDetAMResponse() function in 
DetResponse.c under the LAL package Ref. [1]. 

Table [1] contains the parameters of the injected waveform. On Fig [1] the two polarizations of 
this waveform h^ and hx are shown. The parameters characterizing the source and the detector 
orientation, which are necessary for computing the antenna functions -F+ and Fx are given in 
Table [2j The waveforms appearing at the Hanford and Livingston detectors are plotted on Fig 
[2] (top line), while the bottom line on Fig [2] shows the signals immersed in S5 LIGO-like noise. 
Some characteristic unfiltered and filtered noises can be seen on FigO 

First we calculated the four type of overlaps between the noisy injection and the injected signal 



itself, finding O^ [hn^i,h. 



injection] 



0.8392, 0^[hn,i,h, 



injection 



] = 0.9064, 0^^[hn,i,K 



injection] 



0.8662 and O^J:'^ [hn,i-,hinjecUcm] = 0.8646. Next we defined the following auxiliary quantity for 
all overlaps 

^ ["-n,*! i^template\ 



a 



1 



(4) 



L^ yiT"n,i) i^injectionl 

vanishing for the injected template. We searched for templates with the various as less then 
0.1. 
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Figure 4. The parameters of the injected signal (red) and the templates having a < 1 
according to O^ (green), O^ (yellow), O^^ (dark blue) and all of these (light blue) and when 
the correlated overlap was also taken into account (black). The big dots among the black 
represent those templates that has the smallest values accordint to the correlated overlap. 

The templates were chosen with parameters in the ranges: masses rrii € 3 -^ 10 M©; 
dimensionless spins Xi ^ 0.7 -^ 1; spin angles coskj and ^j random; distance di, angles Q, ip, 
and ^ are fixed identically to the values given in Tables [1] and [2] for the injection, respectively. 
We also assumed as known the time of signal arrival to the detector. Therefore we varied two 
mass and six spin parameters altogether. 



2. Results 

Our analysis is based on matching more than one million templates. The templates with any of 
the cj^, a^ , a^^ < 0.1 were selected and represented on Fig|3]in the parameter planes (mi, 7712), 
(xi5X2) and {cos Ki,ipi) for both spins. The green, yellow and navy dots (for grayscale see the 
legend) represent templates with required values of a^ , errand a^^, respectively; templates 
with all three values of a below the threshold are plotted with larger turqoise dots. Twelve even 
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Figure 5. Best three templates (black) as they would appear at Hanford (left) and Livingston 
(right), compared to the respective injected signals (red). 

larger black dots show templates with all four as (including o"^^^) below the threshold. The 
three largest of them have the lowest value of o"^^^. The parameters of the best three templates 
are shown in Table [3l and they are plotted on Fig [5] as they would appear at the Hanford and 
Livingston detectors. 

Discussion: The masses are reasonably well recovered, although sligtly overestimated with 
any of the a , a , a . The additional monitoring of the correlated match cr^orr imposes 
however a selection effect which improves the estimation of the masses (top left panel of Fig 
H]) . While the recovery of the spin magnitudes is still problematic (top right panel of Fig U]) , 
the estimation of the spin angles seems slightly improved by the use of of the correlated match 
cr^^j.. (Black dots exhibit a belt-like structure on both bottom panels of FiglH) How relevant 
is this feature statistically is currently under investigation. 
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Table 3. The parameters of the templates represented by the four big black points on Fi| 
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